New limit on the present temporal variation of the fine structure constant 
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The comparison of different atomic transition frequencies over time can be used to determine the 
present value of the temporal derivative of the fine structure constant a in a model-independent 
way without assumptions on constancy or variability of other parameters, allowing tests of the 
consequences of unification theories. We have measured an optical transition frequency at 688 THz 
in 171 Yb + with a cesium atomic clock at two times separated by 2.8 years and find a value for 
the fractional variation of the frequency ratio /vb//cs of ( — 1.2 ± 4.4) ■ 10~ 15 yr -1 , consistent with 
zero. Combined with recently published values for the constancy of other transition frequencies this 
measurement sets an upper limit on the present variability of a at the level of 2.0 • 10 -15 yr -1 (la), 
corresponding so far to the most stringent limit from laboratory experiments. 



PACS numbers: 06.20.Jr, 06.30.Ft, 42. 62. Eh 

Motivation to search for a time- or space-dependence of 
fundamental 'constants' has developed because theories 
attempting to unify the fundamental interactions allow 
or even imply variations of the coupling constants 0. 
According to the inflationary model, dramatic changes 
of particle properties took place in the early evolution of 
the universe and it is conceivable that remnants of these 
changes are still observable in later epochs. Variations of 
the constants have been searched for in various contexts 
00. Sommerfeld's fine structure constant a is presently 
the most important test case. As a dimensionless quan- 
tity, a ~ 1/137 can be determined without reference to 
a specific system of units. This point is important be- 
cause the realization of the units themselves may also be 
affected by variations of the constants. 

The only indication for a possible variation of a so far 
comes from an apparent shift of wavelengths of specific 
absorption lines produced by interstellar clouds in the 
light from distant quasars |4J. These observations seem 
to suggest that about 10 10 yr ago the value of a was 
smaller than today: Aa/a = (— 0.54±0.12)-10~ 5 over the 
redshift range 0.2 < z < 3.7, representing 4.7<r evidence 
for a varying a 5| . Assuming a linear increase in a with 
time, this would correspond to a drift rate dhia/dt = 
(6.4±1.4)-10~ 16 yr -1 Recent evaluations of new data 
on quasar absorption lines are consistent with Act = 
0. A stringent limit on the order of \dhia/dt\ < 1 ■ 
10~ 17 yr -1 for the temporal variability of a has been 
deduced from the analysis of a natural nuclear reactor 
that was active 2- 10 9 years ago in the Oklo region in West 
Africa 0,0. Since this analysis is based on the energy of 
a neutron capture resonance, it is not possible to separate 
the influence of a variable fine structure constant from 
possible variations of the coupling constants of the strong 
interaction and of nuclear masses. In unified theories the 
drift rates of these parameters are correlated and it has 



been predicted that nuclear masses and moments may 
show significantly higher relative drift rates than a 0, 
licj . In this sense the limit on d\na/dt from the Oklo 
reactor has to be regarded as model-dependent [llj |. 

In this letter we report a model-independent limit for 
the present temporal variation of a from comparisons 
of atomic frequency standards. We reach a sensitivity 
for Aa/At with At ~ 3 yr that approaches that of the 
quasar observations if a linear change of a is as- 
sumed over a cosmological time span. Possible system- 
atic errors can usually be identified more easily in lab- 
oratory experiments than in astrophysical or geophys- 
ical observations. Precise frequency comparisons have 
recently been performed in order to constrain temporal 
derivatives of fundamental constants 

[HE! (see also 
and references therein). Since these measurements in- 
volve hypcrfinc frequencies, they are sensitive not only 
to a variation of a, but also to a variation of the nuclear 
magnetic moment. In this case a limit on the variability 
of a could only be obtained under the assumption that 
the nuclear properties stay constant, which is difficult to 
justify. 

To avoid the influence of nuclear structure and to ob- 
tain a direct measure of the drift rate of a one can com- 
pare different optical electronic transition frequencies. So 
far, sufficiently accurate atomic optical frequency stan- 
dards have not been operated continuously over extended 
periods of time, but a number of precise frequency mea- 
surements have been performed with reference to cesium 
clocks. For the analysis, we assume that the electronic 
transition frequency can be expressed as 

f = Ry-C-F(a) (1) 

where Ry = m e e 4 /(8ef ) h 3 ) ~ 3.2898 • 10 15 Hz is the Ry- 
dberg constant in SI units of frequency, appearing here 
as the common atomic scaling factor (m e : electron mass; 
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e: elementary charge; eo : electric constant; h: Planck's 
constant). C is a numerical constant which depends only 
on the quantum numbers characterizing the atomic state 
and is assumed not to vary in time. F(a) is a dimension- 
less function of a that takes into account level shifts due 
to rclativistic effects 0, ^5|. In many-electron atoms 
these effects lead to different sensitivities of the transi- 
tion frequencies to a change of a, with the general trend 
to find a stronger dependence in heavier atoms, as was 
first pointed out in |l6j | . The relative temporal derivative 
of the frequency / can be written as 
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The first term din Ry/dt represents a variation that 
would be common to all measurements of electronic tran- 
sition frequencies: a change of the numerical value of the 
Rydberg constant in SI units, i.e. with respect to the 
133 Cs hyperfine splitting. We distinguish here between 
the numerical value of Ry and the physical parameter 
which determines the atomic energy scale. A change of 
the latter would not be detectable in this type of measure- 
ment because the 133 Cs hyperfine splitting is also propor- 
tional to this quantity. Consequently, the detection of a 
nonzero value of 9 In Ry/dt would have to be interpreted 
as being due to a change of the cesium hyperfine interval, 
and in particular of the 133 Cs nuclear magnetic moment. 
The second term in Eq. 2 is specific to the atomic tran- 
sition under study. The value of the sensitivity factor A 
for small changes of a has been calculated for a number 
of cases by Dzuba et al. 0, 0] . Assuming that pos- 
sible drifts of the atomic frequencies are linear in time 
over the period of interest, the drift rate dlna/dt can be 
obtained from at least two measured drift rates of tran- 
sition frequencies 9 In f /dt if the sensitivity factors A for 
the two transitions are different. A plot of dlnf/dt as a 
function of A should yield a straight line where the slope 
is determined by dlna/dt and a drift of Ry would show 
up in the intercept [lOj . 

Optical transition frequencies can be measured very 
precisely in single trapped and laser-cooled ions be- 
cause here perturbations that may lead to systematic fre- 
quency shifts can be controlled very effectively. These 
experiments are well suited for a laboratory search for 
a variation of a because transitions in heavy ions like 
Yb + or Hg + cover a wide range of A values 0, Il5| . 
At PTB an optical frequency standard based on 171 Yb + 
is investigated fl8t Hi) . A single ytterbium ion is 
trapped in a miniature Paul trap and is laser-cooled to a 
sub-millikelvin temperature by exciting the 2 Si/ 2 (F = 
I) — > 2 P\/ 2 {F = 0) resonance transition at 370 nm. 
Two additional lasers at wavelengths of 935 nm and 
639 nm are used to repump the ion from the metastable 
2 D 3/2 (F = I) and 2 F 7/2 levels that may be populated 
during the cooling phase. The electric quadrupole tran- 
sition 6s 2 S 1/2 {F = 0) -> 5d 2 D 3/2 (F = 2) at 436 nm 
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FIG. 1: a) Simplified level scheme of the 171 Yb+ ion. A 
modulation sideband of the cooling laser is used to repump 
the ion from the F = ground state, b) Excitation spectrum 
of the 2 S 1/2 {F = 0,m F = 0) -» 2 D 3/2 {F = 2,m F = 0) 
transition, obtained with probe pulses of 90 ms duration and 
20 probe cycles for each value of the detuning. 



wavelength (688 THz) with a natural linewidth of 3.1 Hz 
serves as the reference transition. The isotope 171 Yb 
has a nuclear spin quantum number of 1/2 so that a 
(F = 0, rriF = 0) — * (F = 2, mp = 0) component of the 
reference transition with vanishing linear Zccman shift is 
available (Fig. la). The reference transition is probed 
by the frequency-doubled radiation from a diode laser 
emitting at 871 nm. The short-term frequency stability 
of this laser is derived from a high-finesse, temperature- 
stabilized and seismically isolated reference cavity. The 
trapped ion is interrogated using the electron shelving 
scheme El with alternating cooling and probe laser 
pulses. A Fourier-limited linewidth of 10 Hz at 688 THz 
and a resonant excitation probability of about 50% is ob- 
tained (Fig. lb). The probe laser frequency is stabilized 
to the atomic reference frequency in a servo system with 
a time constant of ra 10 s. We have observed a fractional 
instability well below 1 • 10~ 15 for an averaging time of 
1000 s. In order to minimize servo errors due to the drift 
of the probe laser frequency (in the range of 0.2 Hz/s), a 
second-order integrating algorithm is used. 

The frequency of the 688 THz 171 Yb + standard was 
measured relative to PTB's primary cesium fountain 
clock |20j| . A hydrogen maser is referenced to the cesium 
clock and delivers a signal at 100 MHz with an instability 
below 2 • 10~ 15 in 1000 s of averaging. The link between 
the maser signal and the optical frequency is established 
with a femtosecond-laser frequency comb generator |2lj |. 
The result of a first series of measurements starting on 
MJD 51889 was 688 358 979 309 312 Hz with a total la 
measurement uncertainty of 6 Hz |22|. The frequency 
was remeasured around MJD 52947 after improvements 
in the resolution of the ionic resonance and in the fre- 
quency comb generator setup. The new value has a total 
lcr uncertainty of 6.2 Hz (see below) and is insignificantly 
lower by 2.3 Hz (see Fig. 2). From these measurements 
we deduce a value for the fractional variation of /yb //cs 
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of (-1.2 ±4.4) • 10 



-15 



yr 
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FIG. 2: Results of absolute frequency measurements of the 
171 Yb + standard versus the date of the measurement (MJD: 
Modified Julian Date). 



The dominant source of systematic uncertainty in the 
m Yb + 688 THz optical frequency standard is the so- 
called quadrupole shift, i.e., the shift of the atomic tran- 
sition frequency due to the interaction of the electric 
quadrupole moment of the -D3/2 state with the gradient 
of the static electric field at the trap center. Such a field 
gradient may be produced by patch charges on the trap 
electrodes and a sensitivity of the transition frequency to 
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field gradients of up to 6 Hz/(V/mm 2 ) is expected 
A static field gradient also leads to modifications of the 
oscillation frequencies of the ion in the trap and an anal- 
ysis of this effect indicates that patch field gradients in 
our trap are smaller than 0.5 V/mm 2 . Other systematic 
frequency shifts are well below 1 Hz and their contri- 
bution to the uncertainty is negligible at present: The 
static magnetic field in the trap is typically 1 fiT during 
the interrogation of the reference transition, leading to a 
quadratic Zeeman shift of 0.06 Hz. The quadratic Stark 
shift due to blackbody radiation emitted from the trap 
electrodes at room temperature is calculated as —0.4 Hz. 
In the absence of a full evaluation, we use an estimate 
of 3.2 Hz for the systematic uncertainty of the 171 Yb + 
frequency standard, dominated by the contribution from 
the quadrupole shift. Frequency comparisons between 
two Yb + ions in two separate traps showed agreement at 
the sub-hertz level [l9|. The total uncertainty of 6.2 Hz 
(corresponding to a relative uncertainty of 9 • 10 -15 ) con- 
tains also a systematic contribution of 2.2 Hz from the 
microwave reference and frequency comb generator and 
a total statistical contribution of 4.8 Hz. 

To obtain the limit on dlna/dt we combine the data 
from 171 Yb + with a published limit on the drift rate 
of a transition frequency in 199 Hg + |l3j. The mercury 
ion is investigated as an optical frequency standard at 
NIST (Boulder, USA). The transition 5d 10 6s 2 S 1/2 -> 
5d 9 6s 2 2 D 5/2 at 282 nm (1065 THz) serves as the ref- 
erence. A first frequency measurement of this transi- 
tion was published in 2001 with a relative uncertainty of 
9 ■ 10~ 15 |23j. A sequence of measurements over a period 



of two years has resulted in a constraint on the fractional 
variation of /H g //cs at the level of 7 ■ 10~ 15 yr -1 (la) 


The transition frequencies of the ytterbium and of the 
mercury ion are the two most accurately known opti- 
cal frequencies and their sensitivities to changes of a are 
quite different 0: A Yh = 0.88 [11] and A Kg = -3.19 
|14j . Since the uncertainties of these numbers are not 
given explicitly, we treat them as being exact, noting 
that changes within 10% would modify our results only 
marginally. Combining the two measurements in a search 
for a temporal variability of a is well justified: The 
Hg + observations cover the time between MJD 51770 
and MJD 52580 so that there is a significant temporal 
overlap. The cesium fountains of NIST and of PTB have 
been compared repeatedly during the relevant period and 
found to be in agreement to within their evaluated un- 
certainties of 1.0 ■ 10" 15 m. 

A further precision test for the constancy of an optical 
transition frequency was recently presented by Fischer 
et al. who measured the IS" — > 2S transition frequency 
in atomic hydrogen in two periods around MJD 51360 
and MJD 52680, using a transportable cesium fountain 
frequency standard |25j . The limit on the fractional drift 
rate for the hydrogen frequency that is deduced from 
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these experiments is (— 3.2 ± 6.3) 
Ah being zero. 

The three limits on the frequency drift rates are con- 
sistent with a linear dependence on A, as assumed in Eq. 
2 (see Fig. 3a). Through a weighted linear regression we 
obtain from the slope of the fit: 
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for the present value of the temporal derivative of the 
fine structure constant at the confidence level of la. Us- 
ing only the data points from the single-ion experiments 
with Hg + and Yb + one finds only a marginal change of 



the result: d\na/dt = (-0.2 ± 2.0) • 10" 



yr 



This 



limit is the result of a model-independent analysis that 
is not based on assumptions about correlations between 
drift rates of different quantities and does not postulate 
the drift rate of any quantity to be zero. To our knowl- 
edge, this is so far the most stringent limit on dhia/dt 
obtained from laboratory experiments. Another recent 
related analysis based on the data from hydrogen and 
mercury is presented in |25l |. 

The intercept of the fitted straight line with the line 
A = in Fig. 3a can be used to determine the drift rate 
of the numerical value of the Rydberg frequency Ry. This 
quantity is of great metrological importance. A nonzero 
value would imply among other consequences that a time 
scale generated with a cesium clock would not be uniform 
with respect to other atomic oscillations. The result that 
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FIG. 3: a) Relative frequency drift rate versus sensitivity fac- 
tor A for the S — > D reference transitions in the Hg + and 
Yb + ions and for the IS — * 2S transition in atomic hydrogen. 
The solid line is the result of a weighted linear regression. A 
significant deviation of the slope from zero would indicate a 
nonzero time derivative of the fine structure constant a. b) 
Constraints on the variability of a and Ry as obtained from 
the experimental points in (a). Stripes mark the la uncer- 
tainty of the individual measurements and the central hatched 
area is bounded by the standard uncertainty ellipse resulting 
from the combination of all three experiments. 



we obtain here is consistent with zero: 
d In Ry 



dt 



-1.6 ±3.2) • 10" 
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The result for the fit based only on the single-ion data is 
dlnRy/dt = (-1.0 ±3.7) • 10" 15 yr -1 . 

Fig. 3b) shows how the measurements of the three 
transition frequencies contribute to the constraints on the 
temporal derivatives of a and Ry. The central hatched 
area is bounded by the standard uncertainty ellipse as 
defined by 
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where i designates the elements (H, Yb, Hg), dYa.fi/dt 
is the central value of the observed drift rate, Ui its lcr 
uncertainty and Xmin = 0-092 the minimized y 2 of the 
fit. 

The kind of analysis that we have applied here may 
be extended to other atomic systems as soon as preci- 
sion data on frequency drift rates become available. Fur- 
ther improvements can be expected since trapped ions 
have the potential to reach relative uncertainties in the 
range of 10 -18 for the comparison of optical transition 



frequencies [T^. The combination of results from differ- 
ent atomic and molecular frequency standards will make 
it possible to disentangle the contributions from the fun- 
damental interactions and thus lead to sensitive tests of 
the consequences of unification theories. 
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